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Abstract 

In this paper, we are concerned with the local existence of strong solutions to the k — e 
model equations for turbulent flows in a bounded domain He R 3 . We prove the existence 
of unique local strong solutions under the assumption that turbulent kinetic energy and the 
initial density both have lower bounds away from zero. 
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1 Introduction 

Turbulence is a natural phenomenon which occurs inevitably when the Reynolds number of flows 
becomes high enough(10 6 or more). In this paper, we consider the k — s model equations [Tj. 1TB] 
for turbulent flows in a bounded domain f2c R 3 with smooth boundary, 


with 


Pt + V ■ (pu) = 0, 

2 

(ри) t + V • (pu <8 u) — A u — V(V ■ u) + Vp = — -V(pfe), 

o 

(ph)t + V • (puh) — Ah = p t + u ■ Vp + Sk, 

(рк) t + V • (puk ) — A k = G — pe , 

(pe)t + V • (pus) - A e = — - 

(p,u,h,k,e)(x, 0) = (p 0 (x),uo(x),ho(x),ko(x),s 0 (x)), 
r)k r)s \ 

Un = (0,0,0,0), 

on on J 


(1.1) 

( 1 . 2 ) 

(1.3) 

(1.4) 

(1.5) 

( 1 . 6 ) 

(1.7) 


/ du l dui\ 2 dir 

^\dxj dxi) 3 ‘ J dxk 


du 1 

dxj 


Pt dp dp 
p 2 dxj dxj' 


( 1 . 8 ) 
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dxj 


dvP_\ __ 2 

ArJ 3 vy 


pA; + m< 


9u fc 


P = P 7 , 


(1.9) 

( 1 . 10 ) 


where Sij = 0 if * ^ j, <% = 1 if i = j, and p, p t , p e , Ci and C 2 are five positive constants 
satisfying p + p t = p e , and it is the unit outward normal to 90. 

The equations (11 . ll) - (ll. 101) are derived from combining the effect of turbulence on time- 
averaged Navier-Stokes equations with the k — £ model equations. The unknown functions p, u, 
h, k and e denote the density, velocity, total enthalpy, turbulent kinetic energy and the rate of 
viscous dissipation of turbulent flows, respectively. The expression of the pressure p has been 
simplified here, which indeed has no bad effect on our study. 

In partial differential equations, k — e equations belong to the compressible ones. In this 
regard, we will refer to the classical compressible Navier-Stokes equations and compressible MHD 
equations, which are also research mainstreams, to carry out our study. 

For compressible isentropic Navier-Stokes equations, the first question provoking our interest 
is the existence of the weak solutions. P. L. Lions mug proved the global existence of weak 
solutions under the condition that 7 > where 7 is the same as in (11.101) and n is the 

dimension of space. Later, Feireisl 00 improved his result to 7 > A) . The condition satisfied by 
7 is to prove the existence of renormalized solutions, which was introduced by DiPerna and Lions 
[6]. When the initial data are general small perturbations of non-vacuum resting state, Hoff jT0| 
proved the global existence of weak solutions provided 7 > 1. The existence of strong solutions 
is another problem provoking our interest in the research of Navier-Stokes equations. It has been 
proved that the density will be away from vacuum at least in a small time provided the initial 
density is positive. If the initial data have better regularity, the compressible isentropic Navier- 
Stokes equations will admit unique local strong solutions under various boundary conditions 
00111 ng. However, when initial vacuum is allowed, it was shown recently in 0 that the 
isentropic ones will have local strong solutions in the case that some compatibility conditions are 
satisfied initially. H. J. Choe and H. S. Kim 0 obtained the unique local strong solutions for full 
compressible polytropic Navier-Stokes equations under the similar condition as in [2]. In 0, the 
technic the authors used is mainly the standard iteration argument and the key point of their 
success is the estimate for the L 2 norm of the gradient of pressure. In the process of studying the 
condition of local solutions becoming global ones, Z. P. Xin (20j proved that the smooth solutions 
will blow up in finite time when initial vacuum is allowed. 

As for compressible MHD equations, the research directions, which mainly contain first the 
existence of weak and strong solutions and second the condition of weak solutions becoming 
strong or even classical ones and the local becoming global ones, are similar to that of Navier- 
Stokes equations. For example, Hu and Wang mmm obtained the local existence of weak 
solutions to the compressible isentropic MHD equations. Rozanova m proved the local existence 
of classical solutions to the compressible barotropic MHD equations provided both the mass and 
energy are finite. J. S. Fan and W. H. Yu in 0 proved the existence and uniqueness of strong 
solutions to the full compressible MHD equations. The method used by J. S. Fan and W. H. 
Yu 0 is similar to that in 0, for example, they are both dependent on the standard iteration 
argument and the estimate for the L 2 norm of the gradient of pressure. 

Under the the hypothesis of the existence of local-in-time smooth solution, the authors of 
[I] prove the existence of small data smooth solution in R 3 . In this paper, we consider the 
local-in-time existence of strong solutions to the k — e model equations (ll.lll - dl.lOl) in a bounded 
domain S2 C R 3 . Our method is similar to that in 0 and 0. However, in the process of applying 
the method to k — e model equations, we find that the regularity of the solutions should be 
higher, which is induced by higher nonlinearity than compressible Navier-Stokes equations and 
compressible MHD equations, than that in 0 and 0. In fact, when we make the difference of 
the n — th and the (n + 1) — th of equation (12.41) and integrating the result, it inevitably comes 
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out the term f djp n+1 djp n+1 ■ /i" +1 . Therefore, we have to use integration by parts, which leads 
to two terms as f ~p n+1 djdjp n+1 ■ h’ l+1 and J p n+1 djp n+1 ■ djh n+1 . Then, by Holder and Young’s 
inequalities, it turns out that ||V 2 p n+1 ||£3 and ||Vp™ +1 \\l°° should be bounded. Thus, we need 
||p||#3 be bounded for a priori estimates. Therefore, from the mass equation enough regularity 
of the velocity field should be imposed. Moreover, due to the strong-coupling property of k — e 
equations, we need corresponding high regularity of unknown functions k and e. 

In a word, the high nonlinearity of k — e equations leads to the necessity of high regularity of 
some unknown functions and thus leads to much difficulties for the a priori estimates. Besides, 
physically, when the turbulent kinetic energy k vanishes, the turbulence will disappear and the 
k — e model equations will degenerate into the Navier-Stokes equations, therefore, without loss 
of generality, we assume throughout this paper that the turbulent kinetic energy k has a positive 
lower bound away from zero , namely, 0 < m < k with m a constant. 

To conclude this introduction, we give the outline of the rest of this paper: In section 2, we 
consider a linearized problem of the k — e equations and derive some local-in-time estimates for 
the solutions of the linearized problem. In section 3, we prove the existence theorem of the local 
strong solutions of the original nonlinear problem. 


2 A priori estimates for a linearized problem 

Using density equation m, we could change (U.lD - dl.lOH into the following equivalent form : 


Pt + V ■ (pu) = 0, 

put + pu • Vu — Ait — Vdivit + Vp = — |V (pk ), 
pht + pu ■ Vh — Ah = pt + u ■ Vp + S k , 
pkt + pu ■ S7k — A k = G — pe, 
pe t + pu • Ve - Ae = 

(p,u,h,k,e)(x,0) = (po{x),u 0 (x),ho{x),ko(x),£ 0 (x)), 

S u -^’ h ’wk>ww) l sn = (0-0,0,0). 

Then, we consider the following linearized problem of m : 


with 


Pt + V • ( pv ) = 0, 

2 

pu t + pv • S7u — A u — Vdi vu + Vp = ——V(p7r), 

o 

ph t + pv ■ V/i — Ah = pt + u ■ Vp + S k , 
pkt + pv ■ Vfc — A k = G — p9, 


pe t + pv ■ Ve — Ae = 


CiG9 

7r 


C 2 p0 2 

7T 


(p, V, h, 7T, 9)(x, 0) = {p 0 (x), U 0 (x), h Q (x), k 0 (x),e 0 (x)), 

(V^,/i,jL, J^)|an = (0,0,0,0). 

\ an an) 


S k = 


G = 


( dv l 

dvi \ 

2 dv k ' 


+ dxi) 



dv l 

8Xn 


+ 


dv l 

dx-j 


/ dv l dv J 
\ dxj dxi 


- -Sij 

3 J 


pn + p, e 


p t dp dp 
p 2 dxj dxj ’ 
dv k \ 1 
dxk) 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 
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where v, n and 9 are known quantities on (0,Ti) x Q with T\ > 0. 

Here we also impose the following regularity conditions on the initial data: 


'0 < m < p 0 , po G H 3 (fl), 
uo g H 3 (Tl), 

< (ho, ko,£o) g H' 2 (n), 



< m < k 0 . 


( 0 , 0 , 0 , 0 ), 


(2.9) 


For the known quantities v, ir, 9 , we assume that u(0) = u o, 7r(0) = ko, 9(0) = Eo and 

suPo^t^dHIff 1 + IMIff 1 + ll^llffO 
+ / 0 T2 (ll 7r llff3 + \\v t \\ 2 H i + ht\\ 2 H i + \\0t\\ 2 H i)dt < ci, 

sup 0 < t <T 2 IMIff 2 < c 2 , 

< sup 0 < t < T2 IMIff 3 < c 3 , (2.10) 

Jo 2 IMI^di < d, 

SUPo<t<T 2 IKIlff 2 < c 5, 

. su Po <t<T 2 IMIff 2 < c 6 

for some fixed constants c* satisfying 1 < Co < Ci(i = 1,2, • ■ ■ ,6) and some time T 2 > 0. Where 

co = 2 + 11(po, «o)||ff 3 + \\(h 0 , ko,£o)\\H 2 - 

And for simplicity, we set another small time T as T=min{co^ 67_16 cj" 10 c^" 8 c^" 8 C 4 2 c$ 2 Cg 4 , T), T 2 } 
and all of the T in section 2 are defined as this. 

Remark 2 . 1 . Here it should be emphasized that throughout this paper, C denotes a generic pos¬ 
itive constant which is only dependent on m, 7 and |f 2 |, but independent of c^ (i = 0 , 1 , 2 , • • • , 6 ). 

Remark 2.2. From the physical viewpoint, we assume that the turbulent kinetic energy k has a 
positive lower bound away from zero , namely, 0 < m < k with m a constant. We do not know 
whether 0 < m < k holds afterwards if its initial value ko > m. 

Next, we would like to prove the following local existence theorem of the linearized system 

m-m- 

Theorem 2.1. There exists a unique strong solution ( p , u , h, k , e) to the linearized problem h2.2A) - 
m> and \2.9\) in [0, T] satisfying the estimates \2.99\) and 12.10(A) as well as the regularity 

p G <7(0, T; H 3 ), p t G <7(0, T; H 1 ), u G <7(0, T; H 3 ) n L 2 (0, T; H 4 ), 
u t G i 2 (0, T; H 1 ), k G <7(0, T; H 2 ) n i 2 (0, T; H 3 ), k t G L 2 ( 0, T; H 1 ), 
e G <7(0, T; H 2 ), E t G i 2 (0, T-,H 4 ),h£ <7(0, T; H 2 ),h t G L 2 (0, T; H 1 ), 

(s/put, y/pkt, \fp£t,\fpht) G L°°(0, T; L 2 ). 

In the following part, we decompose the proof of Theorem 12.11 into some lemmas. 

Lemma 2.1. There exists a unique strong solution p to the linear transport problem f£.2f) and 
( QO|) such that 

vn 

P>—, IMIff 3 (fi) < Cc 0 , Uptll^MO) < <7 c 0 c 2 (2-11) 

for 0 < t <T. 
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Proof. First, applying the particle trajectory method to equation (12.31) . we easily deduce 
p > Po exp ^ - J || Vu|| L =odt^ > po exp(-c 3 T) > y ^ 


and thus 


- < — < C. 

p m 

Second, by simple calculation, we have 

~j^l|p||ff 3 < C||u||h 3||p|| ff 3 +C||V 4 w|| i 2, 
applying Gronwall and Holder’s inequalities, one gets 


I m < 


exp (c J IHIff3dt^ (\\po\\h3+cJ \\v\\ H 4d?j <Cc 0 


for 0 < t < T. 

Next, from equation (12.21) . one obtains 


llPtllff 1 = || V • {pv)\\ H l < C'||/9|| ff 3||u|| i? 2 < Ccoc 2 


for 0 < t < T. 

Thus, we complete the proof of Lemma T2. II □ 

Next, we estimate the velocity field u. 

Lemma 2.2. There exists a unique strong solution u to the initial boundary value problem \2.S\) 
and such that 


WVputWh + M 2 rn + f IIVutHi ids < Cc 5 0 +2 \ \\u \\ H 2 < Ccl +3 \l (2.12) 
Jo 

\\u\\m < Ccf +3 \ic 2 c 5 , ! \\u\\ 2 Hi ds<Ccl +6l c\cl (2.13) 

Jo 

for 0 < t <T. 

Proof. We only need to prove the estimates. Differentiating equation (12.31) with respect to t, 
then multiplying both sides of the result by iq and integrating over D, we derive that 


1 d_ 
2d t 


pu 2 dx + ||Vu t ||| 2 + ||divu t ||| 2 


= - J p t v • V« • u t - J pv t -Vu-Ut-2 J pv ■ Vw f • u t - 
= h + I2 + ^3 + 14 + I5 , 


J Vp t • u t - ^ fWp,)]t ■ u t 

(2.14) 


where we have used equation (12.21) and integration by parts. We will estimate It (i = 1, 2, • • • ,5) 
item by item. 

First, because p has lower bound away from zero, we easily deduce ||u(||l 2 < C\\^fput\\L 2 - 
Therefore, using Holder, Sobolev and Young’s inequalities and (12.101) . we have 


h < C\\v\\ L -\\p t \\ L 3\\X7u\\ L 4u t \\ L e < C\\v\\ L ^\\p t \\ L 3\\Vu\\ L 2(\\y/put\\L^ + ||Vu t || L 2 ) 

< Cc 2 0 4\\Vu\\ 2 l2 + CWy/putWb + (2.15) 
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h < C\\p\\U\\v\\ L ~\\X7u t \\ L 2\\^pu t \\ L * < Ccocl\\^pu t \\ 2 L2 + -||Vut||i a , 


(2.16) 


h <qbllE-lKlUe||v«|| L 3||^u t ||^ <C7 ? - i co||v u ||i3 + r?||^||^ 1 ||Vpw t |li 2 , ( 2 . 17 ) 

where 77 > 0 is a small number to be determined later. 

Next, to evaluate ||Vu ||£ 3 in (12.171) . we can first use Sobolev’s interpolation inequality to get 

||Vu ||£ 3 < C'||V«|| ia ||V«|| L e < C||V U || ia ||V«|| H i. (2.18) 

Then, applying the standard elliptic regularity result to equation (12.31) and using (12.181) . we have 

l|Vu||iji < C^y/putU* + ||«|Ue||V«||i 2 ||V«||J . 1 + ||Vp || i2 + \\S7p\\Mn + IIVttM, 
thus Young’s inequality and (12.101) yield 

IIVuIIhi < Cel 1 {\\yfpu t \\ L 2 + cf||Vu || L 2 + Coci). (2.19) 

Combining (12.171) . (12.181) and (12.191) and using Young’s inequality, we get 

h < C'?? _1 Co 7+1 (ll\/P u tlll 2 + c iII Vu IIl 2 + cgc?) +v\\vt\\ 2 H i\\y/pUt\\ 2 L3- ( 2 - 20 ) 

By integration by parts, we have 

h = JPtdivut < C'c2 _1 ||p t ||z, 2 ||Vu t || L 2 < Cc^c\ + ~||Vu t ||| 2 , (2.21) 


/s 3 


- J p t 7rV ■ u t - (7 Jn t X7p-Ut J pVn t -u t 

< C||p t || i 3||^|| i e||Vn t || i 2 + c4||Vp|| i a||7r t || i a||Vpn t |lL 2 + C'4||V7r t || i 2|| v ^ Mt || i <2.22) 

< Cc$c$c$ + Ctj~ 1 4 + Cr)\\n t \\ 2 Hl \\y/pu t \\ 2 L a + ^||Vu t ||| 2 . 

On the other hand, we easily have 

A J |Vu| 2 = 2 J Vu-Vrq<i||Vn t ||£ 2 + C||Vu||| 2 , 

and 

J \u\ 2 <Cc§\\^u t \\ L 2 \\u \\ L 2 <Cc 0 \\^pu t \\ 2 L 2 +C\\u\\ 2 L 2 . 

Combining (12.141) - (12.161) and (12.20I) - (12.24L we get 
^ t (\\Vpnt\\ 2 L2 + \\u\\ 2 m ) + \\Wu t \\ 2 L2 

< C{ c q 4 + P^cl^cl + vht\\ 2 H i + v\\vt\\ 2 H i)(\\^put \\ 2 L 2 + ||u||) 

+ C(4 7 c 2 c 2 + t I~ 1 Cq' ] ' +3 cI), 

setting i) = and using Gronwall’s inequality, we derive 


(2.23) 


(2.24) 


(2.25) 


II H 1 . 


\WP u t\\ 2 L 2 + IMIffi + [ \\Vu t \\ 2 L2 ds < Cc Q +27 
Jo 


(2.26) 
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for 0 < t < T, where we have used the fact that lim t _>o(|| + IMIffO + Ccg +27 . 

Next, by (12.191) and (12.261) . we deduce 

||Vu||tfi < C'c| +37 c?, (2.27) 


which implies (12.121) by (12.261) . 

Next, we will estimate f ( * ||w||^- 4 df. By the standard elliptic regularity result of equation (12.31) . 
we have 

||V 4 u|| i2 < \\pu t \\ H 2 + ||pv ■ Vu \\ H 2 + \\Vp \\ H 2 + |||v(p7r)|| J?2 . (2.28) 

By simple calculation, the first term of the right hand side of (12.281) can be controlled as 

llp^tllff 2 < C(||p«t||£2 + ||p|| ff2 ||u t || ff 2) < Cco\[ut\\H 2 - (2.29) 

In order to estimate ||V 2 « t || i 2, differentiating equation (12.31) with respect to t yields 


A u t + Vdivrq (2.30) 

_ _ 2 

= ptUt + putt + Ptv ■ Vu + pv t ■ Vm + pv ■ Vut + Vp* + g(V/9 t 7T + p t V7T + V pir t + pVn t ), 
applying the standard elliptic regularity result to (12.301) and using (12.26[) . one obtains 

||V 2 u t || L 2 < C'(||p t || L 4||u t || L 4 + WpuuWl 2 + \\pt\\L*\\ v \\L°°\\Vu\\ L 4 + ||p||i,~||w t ||i4 ||Vt/||x,4 

+IMMMU* + MUlPtWm + ML~\\pt\\m + INMIVttIU* 

+ l|Vp|U4|K|| i 4 + ||p|| i «,||V7r t || i 2) 

< C'(||p Utt || i 2 + Cj +37 C?C 2 C 5 + c| +37 c?||u t |Ui + c 0 c 2 ||ut|| ff i + co||7T t || ff i), (2.31) 


therefore, the key point is to estimate ||/7u tt || i2 . Because we have the fact ||pit tt ||£2 < 
C'|| v /joutt||i2, we could first estimate ||||as follows. 

Multiplying both sides of (12.301) by u tt and integrating the result over yield 


P u lA x + ^ll Vu t|||2 + i^||dfvut||| 2 


= — J PtUt * Utt J Ptv ■ Vti • Utt — J pvt ■ Vtt • Utt - J pu ■ Vu t - Utt - J S7pt • u tt 
— 2 J (Tl'Vpt + P(V 7T + 7T t Vp + pV7T t ) • U t t = J\ + J 2 + J 3 + Ja + J 5 + J 6■ (2.32) 


Using Holder, Sobolev and Young’s inequalities and (12.101) and (12.261) . we get 


Jl < C , C^||/0 t ||l,3||u t || L 8|| v /pU t t|| i 2 < C'C^||p t || L 3(|| V /pU t || i 2 + ||Vu t || L 2)|| v /pU tt || L 2 

< Cclcl\\X/ut \\ 2 L 2 + Uc* +27 c 2 + ^||^u«|||2, (2.33) 


Jl < Cc§\\^U t t\\ L 4Pt\\L4v\\L°°\\Vu\\ L e < Ccl +6 ' 1 ci4 + -^\\y/pU tt \\ 2 L 2, 


(2.34) 


j 3 < Ccl\\^utt\\ L 4vt\\L4Vu\\ Le < Ccg +67 c 4 ||ut||^i + ^11 VPUttWb, (2.35) 
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J 4 < C'c 0 2 ||v|| L ^||ypM tt || L 2 ||VM t || L 2 < Cc 0 c 2 2 \\S7u t \\ 2 L 2 + ^Wy/puttWh, 


(2.36) 


J 5 < Ccl||^«||L 2 ||Vp t || i 2 < Cc 2 J +1 cl + ^y/pattfv, (2.37) 

J& < C'cJ||7r|| i0 o|| v /pu tt || i 2||Vpt||i2 +C'c,5|| v / pM tt || L 2||V7r|| i 4||^|| i 4 

+Cc$ \\y/fm tt \\LA\Vp\\L~>\\nt\\L* + Ccfi \\^pUtt\\L^\\^n\\L 2 
< Ccqc\c\ + Ccq 117r^H^i + - \\y/pu t t\\ 2 L 2, (2.38) 

inserting (I2.33I) - (I2.38I) to (12.321) . then integrating the result over (0,t), we derive 

r [ pu 2 tt dxdt+ IIV^lH, < Ccl + ^c\cl (2.39) 

J o Jn 

where we have used equation (12.31) to get lim t _> 0 l|Vu t (f)||| 2 < Cc^ +i . 

So, combining (12.291) . (12.311) and (12.391) . we obtain 

[ \\pu t \\ 2 H . <CcI + ^c\c 2 2 . (2.40) 

Jo 

In the following, we shall estimate the rest terms of the inequality (12.281) . 

For the second term of the inequality (12.281) . direct calculation yields 

IIpv • Vu \\ H 2 < C||p || ff 2 \\v \\ H 2 ||u||ij3 < Cc 0 C2\\u\\h3, (2-41) 

therefore, we have to evaluate ||ti||. In fact, Applying the standard elliptic regularity result to 
equation (12.31) . we obtain 

||V 3 n || L 2 < C(\\pu t \\ H i + ||pv ■ Vu||ffi + ||Vp|| ff i + ||V(HllffO, (2-42) 

we could estimate the right hand side of (12.421) item by item. 

First, from (12.261) . we have ||ii t || l 2 < Ccq +7 , thus 

\\pu t \\ H i < CcollutlUa + ||Vp||L«.||u t || i2 +C'co||Vu t || i 2 < Cc$ +1 + Cco\\Vu t || L 2 . (2.43) 

Second, using Sobolev’s interpolation inequality and Young’s inequality, we get 
II pv ■ Vu|| ff i < C(\\pv ■ Vu \\ L 2 + || \7{pv • Vn)|| L 2 ) 

< C(co\\v\\ L ~\\Vu \\ L 2 + ||Vp|MMM|Vix||L 2 + c 0 ||Vu|| L 2 ||Vu||J 2 ||V 3 u ||J 2 
+ c oIHIl“I|V 2 m|| L 2 ) < CcJ +3l c\c 2 + ^||u||jj3. (2.44) 

Third, due to (12.111) . we easily derive 

IIVpIUi < Cc 2 0 . (2.45) 

Last, by simple calculation, one gets 

||V(p7r)|| ff i < C'HpH j^ 3 ||7 t||#2 < Cc 0 c 5 . (2.46) 
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Combining (12.391) and (I2.42D - (I2.46D . we deduce 


IMI/ra < Ccq +3 \ic 2 c 5 . (2.47) 

Next, by simple calculation, the third and fourth terms on the right hand side of (12.281) can be 
estimated as 


\\Vp\\m < Cc 3 0 , \\X7(pn)\\ H 2 < Cc 0 \\tt\\ h3 . (2.48) 

Combining (12.261) . (12.28|) . (12.401) . (12.411) and (12.471) - (12. 48p . one deduce 



\\u\\ 2 H At < Ccl +6l c\c |, 


for 0 < t < T. 

Thus, we complete the proof of Lemma [2 .21 

In the following part, we estimate the turbulent kinetic energy k. 


(2.49) 

□ 


Lemma 2.3. There exists a unique strong solution k to the initial boundary value problem iTE51) 
and m such that 


WVpktWh + ||fc|| 2 ffl + r \\Vk t \\ 2 L ,ds < Ccq, (2.50) 

Jo 

(2.51) 


||fc||if2 < Cc$ac%, / \\k\\ 2 H3 ds < Cc 7 0 
Jo 

for 0 < t < T. 

Proof. We only need to prove the estimates. Differentiating equation (12.51) with respect to t, 
then multiplying both sides of the result equation by kt and integrating over f2, we get 

~\\Vpkt\\ 2 L 2 + \\Vk t \\ 2 L2 

= - J p t v -Vfc • kt - j pv t -Vfc • k t - 2 J pv ■ Vfct • k t + J G t ■ k t 

n p 6 


- Pt0-k t - p9 t • kt = ^ Ki, 


" " i=1 

we could evaluate K, (i = 1, • • • , 6) as follows. 

First, using similar method of deriving (12.151) . (12.201) . (12.161) . respectively, one has 

Ki < Ccl4\\S7k\\ 2 L2 + C\\yfpk t \\ 2 L 2 + ^||Vfc t ||| 2 , 


(2.52) 


(2.53) 


K 2 < Cp 1 cl 1+1 {\\^pk t \\ 2 L 2 +c^||Vfc||^2 +c 2 0 c{c\) + p\\vtf H i\\^pktf L 2 , 


(2.54) 


K 3 <Cc 0 c 2 \\^pk t \\ 2 L2 + ^\\X7k t \\ 2 L2 . 


(2.55) 
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Next, differentiating G with respect to t and inserting the result thus obtained to K 4 yield 

K 4 <c J \Vv t \\Vv\\k t \ + C J \ P MVv t \\k t \ + C J \p t MVv\\k t \ + C J IpIWIVuIN 

< C'c||| V / pfc t || i 2||Vu t || i 2||Vu|| i - + CCq || 7 t||loo || Vat ||z, 2 lll/pfct || i2 

+q|7r|| L ~||pt||L3||Vu|U2||fc t || Le +C'c 0 l ||ypfc t || L 2||7r t || i a||Vu|| i 3 (2.56) 

< Cp^coclcl + Cclc\clcl + C\\y/pkt\\v + Cri(\\vt\\ 2 H i + || 7 r t H^i)||||^2 + ^||Vfc t ||£ 2 . 

Last, direct calculation leads to 

^5<||Pt|U3||0|U2||fc t || i6 <Cc^ + C||ypfc t |||2 + ^||Vfc t ||| 2 , (2.57) 


Ke < Cc$\\y/pk t \\ L 2 \\ 6 t \\ L 2 < Crj 1 C 0 + T]\\ 9 t\\ 2 L 2 \\y/pk t \\ 2 L 2 . 
On the other hand, we easily get 

< T||vMi> + c||Vk|||„ 


Combining (I2.52I) - (I2.60I) . we obtain 
d 


\h<Cco\\Vpkt\\h + C\\k\\h. 


(2.58) 

(2.59) 

(2.60) 


dt 


(\\vpkt\\h+\\kr H i)+\\vk t \\i 2 


<C(cl 4 + r) 1 cl 1 + 1 c\ + r]\\v t \\l 
+C(»T 1 c^c^+c^c^), 


vht\\l,+r,\\0t\\h)(\\Vpk t \\ 


L 2 


1 m) 


setting r] = c 1 1 and using Gronwall’s inequality, we deduce 

\Wpktf L 2 + \\k\\ 2 Hl + f pJk t f L 2 &s<Cc\ 
J 0 


(2.61) 


(2.62) 


for 0 < t < T, where we have used the fact that lim^odl+ ||fc||^i ) < Cel 

Then, by the standard elliptic regularity result of equation (12.511 and using (12.6211 . we have 

||Vfc||iri < Ccl\\^pkt \\ L 2 + Cc 0 \\v\\ L ~\\Vk \\ L 2 + C\\Vv\\ 2 L 4 , 


and 


+Cco||7r|| L 4||Vu|| L 4 + Cco\\Q\\l* < Ccfidcl, 


|V 2 fc|| ff i < C(\\pk t \\ H , + ||pv ■ Vfc |Ui + ||G || ff i + \\p0\\ H i). 


To evaluate / 0 * ||fc||^ 3 df, we will estimate the right hand side of (12.6411 item by item. 
In fact, we derive by using (12.6211 and (12.6311 that 

\\pkt\\m < C{\\pk t \\ L * + ||V(pfc t )|M < Cel + Cco\\Vk t \\ L 2 , 


(2.63) 

(2.64) 

(2.65) 


II pv • Vfc||tfi < C{\\pv • Vfc || L 2 + ||V(pu • Vfc)|| L2 ) 

< C(c 0 \\v\\ L ~\\Vk \\ L 2 + ||Vp||^||u|| i 0 o||Vfc || L 2 

+co||Vu|| L 4 ||Vfc|| L 4 +coHU=o||V 2 fc|U 2 ) <Cc$ci4, (2.66) 
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IIG'lUi < C(||Vu|| 2 4 + ||Vv • p ■ 11 L 2 + ||Vw • V 2 u|| L 2 + ||V(Vw • p ■ tt)M 

< C(||Vu||l4 +co|MU-||Vi;|| L 2 + ||Vv||^||V 2 w|U* + co|M|l~||V 2 u|| L 2 
+ ||7r||L“>||V / o||L=o||Vti|| L 2 + Co||Vw||||V7r||_£, 2 ) < Cc 0 ciclc 3 c 5 , (2.67) 


and 


llp^lli ? 1 < C\\p\\ H 3 \\e\\ m < Ccoci. 


( 2 . 68 ) 


Therefore, inserting (I2.65D - (I2.68II to (12.641) and integrating the result thus obtained over (0, £), 
one gets 


||fc||^ 3 d£ < Gcq 


(2.69) 


for 0 < t < T. 

Combining (12.621) , (12.631) and (12.691) , we complete the proof of Lemma 12.31 

In the next part, we estimate the viscous dissipation rates of the turbulent flows e. 


□ 


Lemma 2.4. There exists a unique strong solution e to the initial boundary value problem 
and \2.9\) such that 


\Wp £ A\ 2 l 2 + Ikllff 1 + [ l|V£t|||, 2 <is < Ccq, 

Jo 


\\ £ \\h 2 < Ccq c 2 c 2 


(2.70) 

(2.71) 


for 0 <t<T. 

Proof. We only need to prove the estimates. Differentiating equation (12.61) with respect to t, 
then multiplying both sides of the result by e t and integrating over f2, one obtains 


~\\Vp £ t\\l 2 + \\V £ t\\l 2 


J p t v ■ Ve ■ Et 


pv t ■ Ve • e t - 2 


pv ■ Vet • £t 



(2.72) 


We could evaluate E 4 and E 5 in the first place. Because n has upper and lower bound away from 
zero, direct calculation yields 


Ea<C j (\G t 9\ + \G9 t \ + \GdTT t \)\£ t \ 

< C /(IV* • Vw| + |p t 7rVv| + \pTifS7v\ + |p7rVu t |)|0||e t | 

+C J (|Vw| 2 + |/07rVw|)|fl t ||e t | +C J (|V*;| 2 + |p7rVc|)|0||7r t ||e t | 

< Cch\0\\ L ~\\Vv\\ L oo\\Vv t \\ L 2\\^e t \\ L 2+c4ML'™\\Vp£t\\L 2 \\pt\\Le\\Vv\\ L e\\d\\ Le 

+GcJ|| v / pe t || i 2 ||7rt||z, 6 l|V' l ’||L 6 ll^llL 6 + Gc 0 ||7r|| L » ||0 ||l°° ||'Vu t || £ 21| yfpet\\ L * 
+G||Vpet|| L 2||0t|| i 2||Vu||l oo +Gco||7r|| L »||ype t ||L 2 |l^llLH|Vt;||L~ 

+c4llV^t||w|k t || i .||V«||i.||0|| i -+C'4l|7r||L»||v^ t ||L=i|kt||z,a||V«|| i .||6»|| il . (2.73) 

< Crj-'coclc^clcicl + C4,c\c\cl + Cr}{ ||Vu t ||| 2 + ||7r t ||| 6 + HfltllLOH y/p£t\\h + C\\y/pe t \\ 2 L ,, 
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and 


E 5 <C j \p t 6 2 e t \ + C J \00 t pe t \ + C J \ P e 2 7r t e t \ 

< C\\pt\\L>M 2 M\L* + CcS\\Jpe t \\ L >\\0t\\ L >\\0\\L- + Ccl\\^e t \\ L 4n t \\ L 40\\U2.74) 

< Cf] 1 C 0 Cg + <7Cq4 C 2 + C\\y/p£t\\ 2 L 2 + Cp(\\9t\\ 2 L 2 + \ | 7T t 11 £2 ) 11 sf(>£t 11 2 L 2 + - ||Vet||i 2 - 
Next, using an argument similar to that used in deriving (12.531) . (12.5411 . (12.551) . (12.601) and 


(12.591) . respectively, one gets 

Ei < Cc 2 4 ||Ve||£ a + CWyftetfv + ^||V£ t ||| 2 , (2.75) 

E 2 < C'r?’* 1 Co 7+1 (|| x /p£ t ||| 2 +c?||Ve||^ +c^4) +T?||u t |||i|| x /pe t ||| 2 , (2.76) 

E 3 < Ccocl\\^/~pE t \\ 2 L 2 + —||V£t||| 2 , (2.77) 

^||e||| 2 < C'||£||| 2 + CcoW^petWl*, (2.78) 

and finally 

^||V£||| 2 <i||V£ t ||| 2 +C||V£||| 2 . (2.79) 

Combining (12.72l) - (12.79l) . one obtains 

2 ^(ll ^*||! 2 + liv et |ll 2 


< (7(44 + rj 1 4 7+1 4 + ^IktUff 1 + + 7 7ll 7r t||ff 1 )(llv / P e ‘lll 2 + IkllffO 

+c?7~ 1 co44 c 6 c 3 c 5 + <74444. ( 2 - 8 °) 

setting r] = Ci 1 and using Gronwall’s inequality, one obtains 

WVpetWl* + Iklllrx + f ||V£ t ||| 2 ds < Ccg (2.81) 

Jo 

for 0 <t<T, where we have used the fact that lim t _,.o(|| Hz , 2 + Iklllfi) < Cel- 

Next, applying the standard elliptic regularity result to equation (12(61) and using (12.811) . we 
have 

||V£|| ff ! < C{$\\y/pet\\L* + C 0 |MMV£||l 3 + ||V4||a||4|Le + C 0 || Vv\\ L e ||0|| L e ||7r|Ue + C O \\0\\%) 
< <7(444 + CoCi||V£||f 2 1|V£||J e ), (2-82) 

therefore, by Young’s inequality and (12.811) . one deduce 

INI# 2 < <7444- 

Thus, we complete the proof of Lemma T2.41 □ 
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Finally, we estimate the total enthalpy h. 

Lemma 2.5. There exists a unique strong solution h to the initial boundary value problem \2.4\ ) 
and m such that 


II V~pht\\ 2 L 2 + \\h\\ 2 H i + 


\\Vh t \\ 2 L2 ds<Cc 5 0 , 


I H 2 7 


< Cc, 


1 + 7 . 


,2 2 


Cic: 


(2.83) 

(2.84) 


for 0 < t < T. 

Proof. We only need to prove the estimates. Differentiating equation (12.41) with respect to t, 
multiplying both sides of the result equation by ht and integrating over S 2 , one obtains 


^(11 VphtWb 


I l/i) + II^Mz , 2 


= — J p t v -Vh-h t - J pvt ■ S7h ■ ht — 2 J pv ■ Vht ■ h t + 

7 

+ J u t ■ Vp • ht + J u - X7p t ■ h t + J S kt ■ ht H i- 


Ptt ■ h t 


i=1 


(2.85) 


Firstly, using similar method of deriving the estimates (12.151) . (12.201) and (12.161) . respectively, one 
has 


Hi < c44\\vh\\ 2 La + cwjphtWb + ^\\vh t \\ L *, 


( 2 . 86 ) 


H 2 < Cr) 1 


27+1 

c 0 


( c o c 2 + llv^tlll 2 + c ill^lll 2 ) + r ?ll w t|ll/i II ifpht 111, 2 5 


(2.87) 


H 3 < CcoclW^/phtWl, + ^|| Vh t \\ L *. ( 2 . 88 ) 

Secondly, differentiating equation (12.21) with respect to t yields 

Ptt = ~PtV ■ v + pV ■ v t + v t -Vp + v ■ Vp t . (2.89) 

Therefore, by direct calculation and using (12.89)1 . we derive 

Ha = J [ 7(7 - l )p J ~ 2 Pt - 7P 7-1 G°tV • v + pS7 ■ v t + v t ■ Vp + v ■ Vp t )] • h t 

< Cc 1 - 1 1| p t ||| 4 1| y/pht || i2 + Ccl~ * || pt\\ L s ||Vu|||| y/pht || L ,+C<%-* || JphtW* II 'Vvt || 

II y/~P^t || L 2 ll^t ||l 6 II VpH/,3 -}- Ccq~^ || Vpt ||l 2 ||'1 ? ||l° o || \J~pht || L 2 

< C(c^ +1 ct + r,- 1 # + WVphtWl 2 + v\\vt\\m\\VpM\h) + ^l|V^||| 2 . (2.90) 

Thirdly, simple calculation and (12.261) lead to 

H 5 < C C Q~* \\y/ph t \\ L 2 \\Vp\\ L 3\\u t \\ L 6 < C'Co _ "|| V / p/l t || i 2||Vp|| i 3(||u t || L 2 + ||Vu t || L 2 ) 

< Cco 7 + 1 ||VpM ?,2 + Cco 7+5 + C||Vu t ||| 2 . (2.91) 
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Next, by direct calculation, we know that Vpt = 7(7 — 1 )p 1 2 p t X7 p + 7 p 1 1 V pt- Therefore, 


H 6 <Ccr 2 j \p tM Vp\\h t \ + Ccrf \u\\S7 p t \\h t \ 

< C(%- 2 \\Vp\\ L ~\\p t \\ L 3 \\u\\ L 2 (\\^pht \\ L 2 + hv/^iuo 
+Ccr 1 \\u\\ L A\Vpt\\Li{\\Vpht\\L* + ||Vh t || i2 ) 

< Ccl + 2 ^c 2 2 + C\\yfph t \\ 2 L 2 + ^HVh t ||| 2 . (2.92) 

Last, simple calculation yields \S kt \ < C\ Vu||Vut| + Cp 1 ~ 1 \pt\\S7p \ 2 + C'p 7_1 |V/?t||V/3|, thus 

H 7 <C /|Vr,||Vu| N + Ccr7|p t ||Vp| 2 |/ lt | + Ccr7 |Vp t ||VHN 

< Cc$\\Vv\\ L ~\\Vvt\\ L 2\\y/pht\\ L 2 + Ccq~* \\p t \\ Le WV p\\ 2 L 6 \\y/pht\\ L 2 

+Cco~^\\V p\\ L <~\\V p t \\ L 2 \\^ph t \\ L 2 

< C^coc 2 + c^c 2 + TiWVvtWUy/phtWh + \WphtWl 2 ). (2.93) 

Furthermore, we easily have 

^ t \\h\\ 2 L 2 <Cco\\Vpht\\ 2 L 2 +C\\h\\l 2 , (2.94) 

and 

^||Vh||| 2 < C||Vh||| 2 + ^||Vft,||| 2 . (2.95) 

Consequently, combining (I2.85I) - (I2.95I) . one deduces 

^ t (\\Vpht\\h + \\h\\ 2 Hi) + \\Vh t \\ 2 L 2 

< C(cq 7 +1 C 2 + r/~ 1 Cg 7+1 cf + p\\vt\\ 2 H i)[\\yfpht \\ 2 L 2 + \\h\\ 2 m ) 

+C(c 0 +2 ' 1 4 + r? _1 Co +27 c^c^), (2.96) 

setting p = Cy 1 and using Gronwall’s inequality, we get 

II Vphtfv + \\hf Hl + f ||V/i*||| 2 ds < Cc 5 0 (2.97) 

JO 

for 0 <t<T, where we have used the fact that limt_^o (11 11 ^2 + ||^||^i ) < Cc 5 0 . 

Next, using (12.9711 and the standard elliptic regularity result of equation (12. 41) . one obtains 

||V/i|| ffl < C{4\\Vpht\\ L > + co||w|| L 6||V/i||£3 + c2r 1 ||p t || i2 +^- 1 || W || L a||Vp|| L 3 
+ ||Vu||| 4 + C 7 - 1 ||Vp||i 4 )<Cc| +7 c 2 +Cc 0 c 1 ||Vh||| 2 ||V/ i .||| 1 , (2.98) 

then, Young’s inequality and (12.9711 yields 

\\h\\m < Ccl + \\cl- 

Thus, we have finished the proof of Lemma 12.51 □ 
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Next, let us define Cj (* = !,■■■ , 6 ) as follows: 


ci = C'cq +2t ,c 2 = Cc^ + 3 y c 2 ,c 5 = Ccgaclce = Ccgcfc 2 ,c 3 = Cc 0 ; 


¥+ 3 7„4 


4 9+67 5 2 

C 1 C 2 C 5 ,C4 = Gc 0 'qc. 


-'l u 2> 


then we conclude from Lemma rrn to Lemma [2 .51 that 


{ sup 0 < t < T (lM|tfi + ||fc|| ff i + IklUO 
+ / 0 (Pllff3 + + IIMIffi + lkt||ffi)dt < Ci, 

supo <t<T I \u\\h* < C 2 ,sup 0 < t < T ||u|| ff 3 < c 3 , / 0 T ||u||^ 4 dt < c 4 , 
supo <t<T INU 2 < c 5 ,sup 0<t<T ||e||jj 2 < c 6 


(2.99) 


and 


||p||ff3(Q) < Cco, ||pt||i/i(fi) < Cc 0 c 2 

WVpbtWh + \Mm + Jo WVhtWlids < Ccg, (2.100) 

\\h \\ H 2 < Ccl + 1 c\c\ 


for 0 < t < T. 

Using standard proof as that in [3], we can complete the proof of Theorem 12.11 □ 


3 Existence of strong solutions to the k — £ equations 

Theorem 3.1. There exists a small time T* > 0 and a unique strong solution ( p,u,h,k,e ) to 
the initial boundary value problem \1.1\) - IT.1()\) such that 

p G C(0, T *; H 3 ), pt G (7(0, T*-H 1 ), u G (7(0, T*; H 3 ) n L 2 (0, T*;FT 4 ), 
ut G i 2 (0, T*; H 1 ), k G (7(0, T*-H 2 ) n L 2 (0, T*; H 3 ), k t G i 2 (0, T*; FF 1 ), 
e G C(0, T*; FT 2 ), e t G F 2 (0, T*; Ff 1 ), F G (7(0, T*;H 2 ),h t G F 2 (0, T*; if 1 ), 
(y/put,y/pkt,y/pe t ,y/pht) G F°°(0, T*;L 2 ). (3.1) 

Proof. Our proof will be based on the iteration argument and on the results in the last section 
(especially Theorem 12.11) . 

Firstly, using the regularity effect of classical heat equation, we can construct functions (it 0 = 
u°(x,t),k° = k°(x,t),£° = e°(x,t)) satisfying (u°(x, 0), k°(x, 0), £°(x, 0)) = (uq(x), ko(x), £o(x)) 
and 


( sup 0 < t < r (||u°|| ff i + ||fc°||iri + ||e°||ffi) 

I +/ 0 (Il^°ll/i3 + 11^(11^1 + ||fc°||^i + ||e°||^i)dt < ci, 

I SUp 0 < t < T \\u°\\ H 2 < C2, sup 0 < t < r ||u°|| ff 3 < C3, / 0 T ||u°||^ 4 dt < C4, 

Uup 0 <t<T P°llff 2 < C5, sup 0 < t < T ||e °|| ff 2 < C 6 . 

Therefore it follows from Theorem [2T] that there exists a unique strong solution (p 1 , it 1 , h 1 , fc 1 , e 1 ) 
to the linearized problem (I2.2I) - (I2.6I) with v,i r, 0 replaced by u°,k°,£°, respectively, which sat¬ 
isfies the regularity estimates (12.991) and (12.1001) . Similarly, we construct approximate solutions 
( p n ,u n , h n , fc", e n ), inductively, as follows: assuming that it" -1 , fc" -1 , e" -1 have been defined for 
n > 1, let ( p n ,u n ,h n ,k n ,£ n ) be the unique solution to the linearized problem (I2.2I) - (I2.6I) with 
v,n,9 replaced by it" -1 , fc" -1 , e" -1 , respectively. Then it follows from Theorem 12.11 that there 
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exists a constant C > 1 such that 


sup (||pl*3 + |KIM+ sup (\\u n \\ H3 + \\k n \\ H 2 +\\E n \\ H 2 + \\h n \\ H z) 

0 <t<T 0 <t<T 

+ sup di\^<iu 2 + \\ v ^ k \\ l > + \ wf *?\\ l ’ + uv^riM 

0<*<T 

+ AlKII^ + II^IIIp + llfcril^ + ||e"||ffi + IKII+ \\k n \\ 2 HS ) < c (3.2) 

Jo 


for all n > 1. Throughout the proof, we denote by C a generic constant depending only on m, 7 , 
\n\ and Co, but independent of n. Next, we will show that the full sequence ( p n ,u n ,h n ,k n ,e n ) 
converges to a solution to the original nonlinear problem (11 . Ill - (1 1. 101) in the strong sense. 

Define p n+1 = p— 1 - p n , it— 1 = it— 1 - it”, 7T +1 = ft— 1 - h n , F +1 = ft— 1 - k n , 
£— 1 = £— 1 - £ n , p— 1 = p n+1 - p n = (p n+1 )^ - ( p n y. 

Then, by equations (12.21) - 12.61) . we deduce that (p n+1 , it— 1 , ft n+1 , fc" +1 , e n+1 , p” +1 ) satisfy 
the following equations: 


where 


Pt +1 + V 

• ( r +1 

u n 

+ P n u n ) = 

0, 





(3.3) 

p" +1 u— 1 

+r + x 

+p n+1 u n 

■ Vu— 1 

+ p n+1 u n 

• Vu" + p”u” 

■ Vu” 


-Au n+1 - 

- V(V 

• u 1 

1+1 ) + Vp— 1 = - 

-V(p n+1 k n + p' 
) 

^), 


(3.4) 

P n +X +1 

+ p n+1 h? 

+ p n+1 u n 

. vr +1 

+ p n+1 u n 

■ Vft" + p n u n 

■ V/i" 


— Ah" +1 = 

= pt +1 

+ u n+1 • Vp— 1 +u n 

■ Vp— 1 + 



(3.5) 

p n+ X +1 ■ 

+ p n+1 

K 

+p n+1 u n ■ 

vF +1 

+ p n+1 u n - 

Vfc” 

+ p"u n • 

Vfc” 


-Ar +1 = 

z Gn+1 

- ( 

In - ( P H+1 

E n — p n 

E n ~% 




(3.6) 


„n+l-n+l 

P fc t 
-Ae— 1 


1 —ii+1 1 n+1 

+ P e* + P 


■■C x 


f g : 


r n+jg 


u" ■ Ve 

k n-l 


1 t~ 7 ,-.n 1 jn—n 

It ■ VE + P U ■ VE 


-C 2 


n n+l 


(e n ) 2 


p^e"—) 2 ^ 

k n ~ l )’ 


(3.7) 


5fc >n+ i = + d i u j) - 3 djP n+ 1 d 3 p n+1 , (3.8) 


g!+i = ^<[/x e (^< + diu]) - ^S ZJ (p n+1 k n + Cedin'?)}. (3.9) 

To evaluate ||p” +1 || i 2 , multiplying both sides of equation (13.31) by p— 1 and integrating the 
result over D, we get 

^H/5" +1 ii! 2 = - Jv-(r +1 u n + P n u n )-p n+1 

= -J (p n+1 ) 2 v • it” + p—V ■ VP- 1 + p n p n+1 V ■ u n + p-V • Vp n . (3.10) 

Applying integration by parts to the second term of the second equality of (13.101) and using 
Holder, Sobolev and Young’s inequalities yield 

^Iir +1 |li 2 < g(||V7r n || L oo||p” +1 || 2 2 + ||V^|U 2 ||p - 1 || i 2 + |KIUe||Vp"|Ua||p- 1 IU 2 ) 

< g(l+rr 1 )||p- 1 ||^+^l|Vu"||^ 1 , (3.11) 
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where (13.21) has been used and 0 < 77 < 1 is a small constant to be determined later. 

Next, multiplying both sides of (13.41) by u n+1 and integrating the result thus derived over il, 
one obtains 


Id 

2dt 


\\V^u n+1 \\h + \\Vu n+1 \\li + IIV • u n+1 \\b 


-2 


p- ' * u■{ ■ u "' * - / p n+ 1 u n ■ X7u n ■ u n+1 - J p n u n ■ X7u n ■ u n+L - 

5 


■n+1 —n+1 


X7p n+1 • u 


+ / —V(p n+ 1 k n +p n k")-u 


T7 n + 1 — 




(3.12) 


i= 1 


Using Holder, Sobolev and Young’s inequalities and (13.21) . we estimate L\, L 2 and L 3 , respec¬ 
tively, as follows: 


Li < C||p" +1 || i2 |K|| i 3K +1 || i 6 < C||p"+ 1 || L2 |K|| i 3(||y^Tn"+ 1 |U 2 + ||Vu ,l+1 || i 2) 

< 5||<||| 3 ||p” +1 ||i 2 +U||v / ^ TT M” +1 |li 2 + i||VM n+1 ||^ 2 , (3.13) 


L 2 < C\\p n+1 \\ L 2 \\u n \\ Le II Vu n || L e ||m” + 1 || L e 

< C\\r +1 \\l> +C\\Jp^u n+1 \\l. + i||VU" +1 ||| 2 , (3.14) 


L 3 < C\\u n \\ Le \\Vu n \\ L 4^/^ TI u n+1 \\ L l < 57+11 ^/^+ T > +1 ||| 2 +v\\u n \\ 2 Hl . (3.15) 

And then, one deduces by integration by parts that 

p n+ 1 \7-u n+1 <C Jp n+ 1 V-u n+1 < C\\p n+1 \\ 2 L 2 + i||V7l" +1 ||| 2 , (3.16) 



L 5 = IJ p n+ 1 k n X7 ■ u n+1 - k n Vp n ■ u n+1 - p n X7k n ■ u n+1 

< C\\p n+ 1 \\ L 2 \\S/u n+1 \\ L2 + C\\T\\ L e\\\'p n \\ L 3 \\y/^U n+1 \\ L 2 + C\\S/T\\ L 2 \\y/^U n+1 \\ L2 

< C( 1 + r ? - 1 )(||p" +1 ||i 2 + ||v / ^ TT ^ +1 |li 2 ) + ^l|VU n+1 Hi 2 + Cp 11+1 + . (3.17) 


Inserting (13. 13D - 13.171) to (13.121) and using inequality ||M ra+1 || i2 < C\\y/p n+ 1 u n+ 1 \\ L 2 , one has 

^||+^u"+|++K +1 || + (3.18) 

< C( 1 + p- 1 + IKII|3)(||p" +1 ||| 2 + II v / ^ t u” +1 ||| 2 ) + C v \\Tf H1 + Cp\\u n \\ 2 m . 

7~l~ j - X 

Then, multiplying both sides of (13.51) by h and integrating the result thus got over f2, one 
obtains 


^n^^ +1 ii! 2 + iiv7r +1 n| 2 

= - J p n+1 K ■ ^ +1 - J P n+1 u n ■ Vh n ■ h n+1 - J p n u n ■ \7h n ■ h n+1 


(3.19) 


ip\ 


wfl+1 


-u n+1 


V+ +1 + u n 


vr +1 ) 


—n+1 

h + 


(*^7c,n+l *^fc,n) ‘ ^ 


■n+1 


E M - 
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First, using similar method of deriving (13.131) . (I3.1dl) and (13.151) . respectively, one easily obtains 


m x < c\\K\\Up n+1 \\i* +qiv 9 ^r + 1 n 2 L2 + ^nvr + 1 m 2 , (3.20) 

m 2 < c\\r +1 \\h+c\\V^v i+1 \\i 2 + l||vr +1 ni 2 , (3.21) 

M 3 < C77- 1 ||v / ^ T r + 1 ||i 2 + p\\u n f m . (3.22) 

Second, simple calculation leads to 

M 4 = J [ 7 (p” + 1 ) 7 -V +1 - 7(p n ) 7- Vr] • h n+1 + J u n+1 ■ Vp n+ 1 h n+1 
+ J u n ■ Vp n+ 1 h n+1 . (3.23) 

By the differential mean value theorem, the first integral of (13.231) can be controlled as 

J b c^+y-vr 1 - 7(p n ) 7 “V] • h n+1 

— C J \p n+ 1 \\prv +1 \ + j 7 (p") 7 " 1 ^ 1 • r +1 . (3.24) 


By equation (13.31) . the second integral on the right hand side of (13.241) can be estimated as 

J 7 (p n ) 7 ~ 1 pr +1 • h n+1 = - j 7 (p”) 7 “ 1 V ■ (p n+1 u n + p n u n ) • h n+1 
<C j \Vp n \\h n+1 \\p n+1 \\u n \ + C J |p" +1 |K||V7T +1 | 

+C J (|Vp n ||U n | + |p n ||Vu n |)|ft n+1 |. (3.25) 

Then, the second integral on the right hand side of (13.231) can be controlled as 

J u n+1 • Vp n+i r +1 <C J |tZ" +1 ||Vp" +1 |r +1 |. (3.26) 

Next, applying integration by parts to the third integral on the right hand side of (|3.23l) . we 
easily get 

J u n ■ Vp n+1 h n+1 <C j\Vu n \\p n+1 \\h n+1 \ +C J K||p" +1 ||Vr +1 |. (3.27) 

Consequently, combining (13.23D - (13.27D and using Holder, Sobolev and Young’s inequalities and 
m , one obtains 

M 4 < C(l + r ? - 1 )(||p" +1 ||i 2 + \\sf^h n+1 \\ 2 L .) 

+Jir +1 n^ + ^iivr +1 ni 2 +c v \\u n r H1 . ( 3 . 2 s) 
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Finally, we evaluate M 5 . Direct calculation yields 


M 5 <C y'(|Vu T, | + |V« n - 1 |)|Vu n ||r +1 |+C J |p n+1 ||Vp n+1 | 2 |*" +1 | 

+/ • ®“ +I + / 

<Cy'(|VM"| + |Vw"- 1 |)|Vu ri ||ft” +1 | + C' J |p" +1 ||Vp n+1 | 2 |r +1 | 

+U J \X7p n \\X7p n+1 \\p n+1 \\h n+1 \ + C j |V 2 p n+1 ||p” +1 ||^" +1 | 

+c J ivp" +i iip" +i nvr +1 i+c' J\vp n \ 2 \\r + 1 \\h n+1 \ 

+C J \V 2 p n \\-p n+1 \\h n+1 \ + C J |Vp”||p" +1 ||V7T +1 |. (3.29) 


Then, applying similar method of deriving (13.281) . one deduces 

M 5 < 5( 1 + rr 1 )(||p" +1 || 2 2 + II +v\\u n \\h + ^||V7T +1 || 2 2 . (3.30) 

Consequently, inserting (13.20l) - (13.22|) . (13.281) and (13.301) to (13.191) . one gets 

A||y^Trr +1 ii 2 2 + iir +1 n 2 fl 

< 5(i+ v- 1 + ii/ l riii3)(iip n+i ii 2 L2 + 11 v / ^ T r +1 ni 2 ) 
+\\\u n+1 fm+Cv\\u n \\ 2 H i. (3.31) 

For the turbulent kinetic energy k , using similar method of deriving (13.191) . one easily deduces 
from equation (13.61) that 


^nVp^ n+1 ni 2 + iivr +1 n| 2 = ~J r +l K ■ k n+1 - j r +1 u n • • k n+1 ( 3 . 32 ) 


r n+1 


r n+1 


- / P n u n ■ S7k n • k" 1 * + / (G n+1 -G n )-k—- I ( P n+1 e n - p n e n ~ 1 ) • k n+ = Y, Ni- 

2 — 1 


We first evaluate N±. Using inserting items technic, one easily gets 
N 4 <cJ (|Vu n | + |Vu n - 1 |)|Vlt Tl ||fc ,l+1 | 

+C J (|Vu"| + |V?i ra_ 1 ||p" +1 | + |V«"- 1 ||r|)|r +1 |. (3.33) 


Using Holder, Sobolev, and Young’s inequalities and (13.21) . we have 

^4<C(l+77 1 )(||p n+1 || 2 2 + || \/p n+1 k || 2 2 ) + Cpllfc Wh 1 + Cp\\u n \\ 2 H i. (3.34) 

Second, we estimate N 5 . Using similar method of deriving (13.331) and (13.34|l . we have 

N 5 =J (p" +1 £ n + p n £ n ) ■ k n+1 < C(||p" +1 || i2 || e "|| L oo + Ir|| i a||p”|| i 3)||y^r +1 || i 2 

<5(1 + 0(11 V^ TT r +1 || 2 2 + ||p" +1 ||| 2 ) + 5r ? r lllri- (3-35) 
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Next, using similar method of deriving the estimates of (13.131) . (13.141) and (13.151) . respectively, 
one easily gets 

Ni < C\\k?\\h\\r +1 \\h +C\\V^ T k n+1 \\ L . + i||vF +1 || 2 i2 , (3.36) 

N 2 < C\\r +1 \\ 2 L*+C\\V^T +1 \\ 2 L2 + i||vF +1 || 2 2 , (3.37) 

N 3 < Cr ? - 1 ||V^r +1 ||| 2 +^||B"||^. (3.38) 

Consequently, inserting (I3.34D - (13.38D to (1 3.32 [) . one deduces 

^iiv / ^ T r +1 ii| 2 + nr +1 iiF (3.39) 

< c( 1 + v- 1 + \\k?\\hm v / F TT r +1 ||| 2 + \\r +1 \\h) + Cr,(\\Tf m + IKIIF + iniF). 

Next, multiplying both sides of m by e n+1 and integrating the result over O, one gets 


Id 
2 dt 


\\VJ^ TI £ n+1 \\ 2 L 2 + \\Ve n+1 \\ 2 L2 = - J p n+1 e? -e n+1 - J -p n+1 u n ■ Ve n • £ n+1 

- / P n u n • V £ n ■ £ n+1 + Cl 


r* c- n n c- n ~ i 
^71+ l fc V 


-n+1 


-c 2 


r p" +1 (e") 2 pFC 1 " 1 ) 2 ! 


fc 71 - 1 


■£ n+1 


— X] Qi' 


(3.40) 


i= 1 


Using an argument similar to that used in deriving (13.131) . (13.141) and (13.151) . respectively, we 
obtain 

1, 


01 < c|| £ ril!3|IF +1 lli 2 + c\\^ TT e n+1 \\l* + ^IIV£" +1 |l! 2 , 


Q2<C\\p n+1 \\h+C\\V^ TT e n+1 \\l* 


|Ve" +1 ||| 2 , 


(3.41) 

(3.42) 


Oa < (3.43) 

Next, direct calculation leads to 

04 < c j (|Vu n ||Vu"| + |Vu"||Vu n - 1 |)|£ n ||£ n+1 | + C J (FI + |£ n - 1 ||F|)|Vu"- 1 | 2 |e n+1 | 

2Ci r f {d j ufp n+1 k n £ n k n ~ 1 - d j u?- 1 p n k n - 1 e n ~ 1 k n ) _„ +1 

“ 3 F/ FF^ £ 

< J (|Vu"||Vu ra | + |Vu"'||Vu n - 1 |)|£ rl ||£ rl+1 | 

+c J (F| + |£"- 1 ||F|)|Vu"- 1 | 2 f +1 | 

+C J (|Vu"| + |Vu n_1 ||F +1 | + |Vu Il - 1 ||F|)|£ n ||£" +1 | 

+C J (F| + |£"- 1 ||F|)|Vu"- 1 |F +1 | (3-44) 

< C(1 + F 1 )(llv / F TT e n+1 lli 2 + IIF +1 lF) 

+Cr?(||u n ||| f i + \\k \\ 2 h1 + IFIlF) + -||V£ rt+1 || 2 2 . 
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Finally, using similar method in deriving the estimate of Q 4, one deduces 


1 , 


Qs < C(1 + 0(11y/OOll* + \\p n+1 \\h) + Cv\\Ve n f L2 + -||Vs" +1 |ll- (3.45) 


Consequently, inserting (I3.41D - (13.45D to (13.401) . one derives 

d_ 
df 1 


-11 Voodoo 1 " 2 


ff 1 


(3.46) 


< C( 1 + ?7 _1 + Ik?||i 3 )(|| n/OOIIs + ||p” +1 ||| 2 ) + C V (\\k n \\ 2 m + K||^ + imO 


In the end, combining (Id.111) . (Id.181) . (Id.dll) . (13.391) and (Id.461) and setting ip n+1 (t) = 

'r-n+l. 

I ^ I II . / nUTi ff. II 

L2 


\r +1 \\h + llx/O^OH. + IIV'OOK + ll\/Ofc n+1 ||| 2 + llvOOO, we get 


(3.47) 


f t P n+1 (t) + \\u n+ Tm + KHlff 1 + lOK + \\e n+1 W 2 m 
< c(i +k + Kills + n^riiia + Kills + kiioo*) 

+Cr](\\u n \\ 2 H i+\\k \\ 2 h i + ||e n |||-i). 

Setting I% (t) = C( 1 + 17“ 1 + Kills + Kills + Kills + Kills) and applying Gronwall’s 
inequality to (13.47P yield 


<p n+i (t)<C V 


ex P / ^(s)ds 


Klliff 


life 


Iff 1 


+ P 


■n||2 

Iff 1 


)ds 


(3.48) 


where it should be noted that cp n+1 (0) = 0. 
Since 


J"(s)ds < Ct + Cr)-H + C, 


(3.49) 


setting T < r] < 1, then we have 


J"(s)ds < CC 


(3.50) 


for t <T. 

By (ld.48D - (13.50D . integrating (13.471) from [0,i], one derives 


0(0 f 

Jo 


POl I 1 


nl-ra+l|,2 

\\h III 1 


+ l|fc” +1 ||^ 


OllO* 


< CCp 


Kllff. 


II k 


Iff 1 


+ e 




)ds 


J”(s)ds exp / I™(s)ds ) + 1 


< Cr]exp(C) 


Kllff 1 


life 


Iff 1 


+ r 


iff 1 


)ds 


(3.51) 


for T* := min{T,T}. 
Therefore, we have 


^2 supo<t<T<p n+1 (t) + ^2 


lOllff 1 


\\r +1 \\ 2 m 


^ +1 \\ 2 m + \r +1 \\ 2 m)ds 


< Cr] exp(C) ^2 / (ll u "llffi + ll fc "llffi + ll £ "llffi)ds. 


(3.52) 
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Thus, choosing 77 such that Cryexp(C) < h, one deduce 

OO OO 

^sup 0 <t< T if n+ 1 (t) + ^2 / P" + 1 |lffi ds 

n =1 n=l 

1 00 „t 

+oE / (K +1 n^ + r n^ + ir +1 ii^)d« 

2 n—1 ^ 

<CC < 00 . (3.53) 

Therefore, we conclude that the full sequence (p ra , u n , h n , k n , e n ) converges to a limit (p, u, h, k, e) 
in the following strong sense: p n —> p in L°°(0,T; L 2 (fl)); (u n ,h n ,k n ,e n ) ( u,h,k,e ) in 

L 2 (0, T; H 1 ^)). It is easy to prove that the limit ( p,u,h,k,e ) is a weak solution to the original 
nonlinear problem. Furthermore, it follows from (13.21) that ( p,u,h,k,e) satisfies the following 
regularity estimates: 

su p (IIpIIh 3 + llptllff 1 ) + su p (IMI-f/ 3 + II^IIff 2 + Iklli? 2 + II^IIf/ 2 ) 

0 <t<T* 0 <t<T* 

+ sup (\\Jpu t \\ L , + \\y/pht\\L* + WVphh* + \Wpet\W) 

0<t<T* 

+ [ T (IMIff* + \\h t \\ 2 m + \\k t f m + \\s t f m + \\u\\ 2 Hi + \\kf H3 ) < c < OO. 

Jo 

This proves the existence of strong solution. Then, we can easily prove the time continuity of the 
solution (p,u,h,k,s) by adapting the arguments in [2j[5]. Finally, we prove the uniqueness. In 
fact, assume that (pi, u±, hi, k\, £ 1 ) and (p 2 , U 2 , ^ 2 , k 2 ,s 2 ) are two strong solutions to the problem 
(II. ID - (II.10H with the regularity (13.IF Let (p,u, h, k,e) = {pi — P 2 , u\ — U 2 , hi — / 12 , k\ — fc 2 , £1 — £ 2 ). 
Then using the same argument as in the derivations of (13.1 IF (13.181) . (I3.31F (13.391) and (I3.46F 
we can prove that 

^ t m\h + + \\Vp^\\h + wvp^wh + wv^nh) 

< R(mm\h + wvp^wi* + wvp^wi* + iiv^iii* + wv^nh) 

for some R(t ) € L 1 (0,T*). Thus, by Gronwall’s inequality, we conclude that ( p,u,h,k,e ) = 

(0, 0,0, 0,0) in (0, T*) x 12. This completes the proof of Theorem 3.1. □ 
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